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Abstract 

A new finite lattice calculation of the low lying bound state energies in the 
massive Schwinger model is presented, using a Hamiltonian lattice formula- 
tion. The results are compared with recent analytic series calculations in the 
low mass limit, and with a new higher order non-relativistic series which we 
calculate for the high mass limit. The results are generally in good agreement 
with these series predictions, and also with recent calculations by light cone 
and related techniques. 

PACS Indices: 12.20-m, 11.15.Ha 
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I. INTRODUCTION 



The Schwinger model [|I],0], or quantum electrodynamics in two space-time dimensions, 
is the simplest of all gauge theories. It has many properties in common with QCD, such as 
confinement 0, chiral symmetry breaking, charge shielding 0, and a topological 6*- vacuum 
0-0. For these reasons it has become a standard test bed for numerical techniques designed 
for the study of QCD, and has been the subject of intensive study over the years in order 
to provide some insight into QCD. 

The purpose in this paper is to carry out a new Hamiltonian finite-lattice calculation of 
the 'positronium' bound-state energies in the massive Schwinger model, which can be com- 
pared with some recent analytic perturbation theory calculations. The massless Schwinger 
model is exactly solvable, and is equivalent to a theory of free, massive bosons A mass 
perturbation theory expansion about the zero-mass limit can be performed to treat the case 
of low-mass fermions 0. This expansion has recently been carried to second order by Vary, 
Fields and Pirner ^ and by Adam Q]. In the case of large fermion mass an expansion 
about the non-relativistic limit can be performed, which involves a Schrodinger equation 
with a linear Coulomb potential . The next-to- leading terms in this expansion have been 
discussed by Coleman and are evaluated in this paper. 

There have been many numerical calculations of the bound-state spectrum. In the Hamil- 
tonian lattice formulation, strong-coupling series were first calculated by Banks, Kogut and 
Susskind extended by Carroll et al. and have recently been taken to tenth order 
in a; = Xjg^a^ by Hamer, Zheng and Oitmaa [|ll|]. The series can be reliably extrapolated 
to the continuum limit, and give detailed and accurate information about the spectrum, 
although not quite as accurate as the finite-lattice method presented in this paper. Finite- 
lattice Hamiltonian calculations were performed by Crewther and Hamer |]12| and Irving and 
Thomas over fifteen years ago. Hamer et al. [|Tl|] did some new finite-lattice calculations, 
but used free boundary conditions; here we use periodic boundary conditions, which should 
give better convergence. 

The model has also been a showcase for the light-front approach to field theory. An 
early and quite accurate variational calculation in the infinite momentum frame was carried 
out by Bergknoff [^. Discrete Light-Cone Quantization (DLCQ) has been applied to the 
Schwinger model by EUer, Pauh and Brodsky [^, and produced good results, not only for 
the lowest states, but for a wide range of excited states. Mo and Perry [|16| have used a light- 
front Tamm-Dancoff approach using up to four-body states which produced outstanding 
results. A calculation using fermion mass perturbation theory with a variational approach 



was carried out by Harada et al. [|T^] to O(m^). Harada et al. ||T8[ showed that the six- 
body contribution to the lowest-lying states was negligible. Kroger and Scheu |TP| use a 
momentum representation with a lattice Hamiltonian corresponding to a 'fast moving frame' 
which also gives good results. Berutto et. al. |]2D| use a improved strong coupling expansion 
which shows a dramatic improvement of physical parameters compared to previous strong 
coupling expansions. 

In Section II of this paper the continuum formalism and the known analytic results for 
the spectrum of the Schwinger model are reviewed. The extension of the non-relativistic 
series to order (g/mY^^ is outlined, and the Hamiltonian lattice formulation is summarized. 
In Section HI the methods of calculation are outlined briefiy, and in Section IV our results 
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are presented. It is our opinion that these are the most accurate numerical estimates of the 
low-lying bound-state energies yet obtained. Our conclusions are summarized in Section V. 



II. THEORY 

A. Continuum Hamiltonian 

The massive Schwinger model is QED in (1+1)D, and can be defined by the Lagrangian 
density 

£ = -^F^.F^'^ + ^ _ ^ 4 _ (1) 
where ip is a two component fermion field, and 

F^. = d^A, - d,A^ (2) 

The coupling g in (H-l) dimensions has the dimensions of mass. Choosing the time-like 
axial gauge: 

Ao = (3) 

the Hamiltonian density is found to be: 

n = -ii)-i\di + igAi)ij + m^ij + ^F^, (4) 

The stress-energy tensor F^" has only one component in one spatial dimension: 

= -A^ =E (5) 

The remaining gauge field component is not an independent degree of freedom, but can 
be eliminated by the equation of motion (Gauss' law): 

diE = -diA^ = gtp-f'^ip (6) 

Equation @ determines the electric field E up to a constant "background" field term P|. 
The background field will be set to zero for this study. 

In the massless case, the theory has been solved by Schwinger and becomes equiv- 
alent to a theory of free, massive bosons, with mass 

Ml 1 

— = ^ ~ 0.564 7 

9 V^r 
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B. Expansions around the massless limit 

For small electron mass m/g, one can obtain analytic estimates by perturbing about 
the massless limit. Caroll, Kogut, Sinclair and Susskind found that the lowest-mass 
( "vector" ) state has mass 

Ml 1 ^, 

— = — + - + ■ ■ ■ ^ 0.564 + 1.78 - + ■ ■ ■ , (8 

9 Vt^ g 9 



while the ratio of the next-lowest ( "scalar" ) mass to the vector mass was 

^ = 2 - 27rV^(-)2 + . . . ~ 2 - 197(-)2 + ■ ■ • , (9) 
Ml g g 

where 7 ^ 0.5772... is Euler's constant. 

These results have been extended to second order by Vary, Fields and Pirner |^ and by 
Adam § 

— = 0.5642 + 1.781 — + 0.1907 — f + ■ ■ ■ 10) 

9 9 9 

Adam p| also found 

- = 2 - —{-f + 2 - 24.625(-)^ + ■ • • (11) 

which differs from the result of Carroll et al. by a factor of 8. Hence 

^ = 1.128 + 3.562(-) - 13.512(-)2 + ■ ■ • (12) 
9 9 9 



C. Non-Relativistic Expansion 

In the non-relativistic limit the bound states can be described by a Schrodinger equation 
with a linear Coulomb potential |P : 

i— + -9'^\x\\mix) = E^ix) (13) 
\m 2 J 

where the non-relativistic energy is, 

E = M -2m (14) 

Using reduced variables 

z = i^f'^ (15) 
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equation ([T3|) becomes 



(^-N + AMz)=0 



(16) 



which is Airy's equation. The solutions are 

where the eigenvalues are given for the symmetric states by 

Al'{-\n) = 

and for anti-symmetric states 

Ai(-A„) = 0. 
Hence the lowest vector (symmetric) state is, 

— ~ 0.642(^)1/=^ as m/g oo 
g m 

and for the lowest scalar (antisymmetric) state 

1A73(—Y''' as m/g ^oo. 



(17) 



(18) 



(19) 



(20) 



(21) 



E2 , 5^x1/3 

g m 



(22) 



Coleman has calculated the next order terms in the large-mass expansion. Using his 
results, the reduced Schrodinger equation becomes 



IT \mj 



+ 



V4m/ 



ij{z) = 



(23) 



The first extra term is merely a constant self-energy term for the fermions, while the effect 
of the remaining terms can be found to leading order by taking their expectation value with 
respect to the leading-order eigenfunctions. By including these terms the results are 



Ei/g = 0.6418 



■^^1/3 _ 



m 



0.25208 

n m "m 



9_-.5/3 



and 



E2/g = 1.4729(^)1/3 - -- + 0.10847(^)^/3 
m n m m 



(24) 



(25) 
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D. Lattice Hamiltonian 



Using equation for the Hamiltonian density, a mapping onto the lattice can be 
achieved. The "staggered" lattice approach of Kogut and Susskind |^ is used. Single 
component fermion fields on site n are defined to obey the following anti-commutation 
rules: 



{0l'(n),0(m)} = (5, 
{0(n),0(m)} =0 



nm 



(26) 



The sites n to n + 1 are connected through a link operator: 

U{n, n + 1) = e^"'^^i(") = e^^^'^) (27) 
The equivalent lattice Hamiltonian to equation (^) is: 

A'' N 2 N 

H = -^J2 [0Hn)e*'("V(^ + 1) - h.c] + m (-l)"0t(^)0(^) + ^ ^ L^{n) (28) 

n=l n=l ^ n=l 

The number of sites N is even and the correspondence between lattice and continuum fields 
is: 



and, 



—9(n) -A\x) 
ag 

(30) 

gL{n) ^ E{x) 
The gamma matrices are represented by, 

A "compact" formulation has been chosen where the gauge field is an angular variable 
on the lattice and L{n) is the conjugate "spin" variable, 

[0{n),L{m)]=t6n^^ (32) 

so that L{n) takes on integer eigenvalues L{n) = 0, ±1, ±2... (which mimics the quantization 
of the flux in one dimension []TU[ in the continuum Schwinger model). Banks et. al. |[T^ 
have demonstrated that the lattice Hamiltonian is gauge invariant and reproduces the Dirac 
equation and QED in the continuum limit. 

The one component fermion fields can be replaced by Pauli spin operators via a Jordan- 



Wigner transformation |TD 
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<P{n) = Y[[ta,il)]^-in) (33) 

l<n 

Hn) = UHcr-M^^n) (34) 

l<n 



giving 



j:WHn)e^'(-^a-in + 1) + h.c] + lm$:(-l)V3(n) + L\n) (35) 

Now define the dimensionless operator, 

W = -^H = Wo + xV (36) 



where 
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V = ^[a+(n)e^^("V-(n + 1) + h.c. 



Wo = Y: L\n) + I E(-l)"^3(r^) + iV/./2 (37) 



For X <^ 1 strong couphng perturbation theory can be used on the modeL In the strong 
couphng hmit the unperturbed ground state |0) has, 

L{n) = 0, asin) = -(-1)", all n (38) 

and the operators exp{±i6{n)) act to raise/lower the electric field L{n) by one unit, according 
to equation (^). Using Gauss' law the gauge field can be ehminated, 

L(n)-L(n-l) = i[a3(n) + (-!)"] (39) 

Using periodic boundary conditions, L{N) = L{0), there remains one independent gauge 
degree of freedom. This is called the "background" electric field which results in half- 
asymptotic particles 0. The ground state energy is the lowest lying eigenvalue, uq, in the 
sector containing the unperturbed ground state |0). The first excited state, c^i, is the lowest 
eigenvalue among the "vector" states, corresponding in the strong coupling limit to 

1 N-l 

\l) = ^T. [fx+(n)e''("V-(n + 1) - h.c] |0) (40) 



The second excited-state energy uj2 is the lowest of a "band" of excited states in the vacuum 
sector, corresponding in the strong-coupling limit to the state. 



1 N-l 

|2) = ^ E [cy^{n)e^'^^^a-{n + 1) + h.c] |0) . (41) 
viV „=i 
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III. METHOD 



The Hamiltonian matrix in the spin representation of equation ( P7D has "^^Cn basis 
states, times the number of background field values allowed. The matrix for each finite 
N was solved by exact diagonalization using the Conjugate Gradient Method (CGM). The 
basis states were 'symmetrized' with respect to translations, assuming periodic boundary 
conditions. The simulations were run for lattices up to = 22 sites using approximately 
1.5 gigabytes of memory and storing approximately two million states. The "vector" and 
"scalar" states can be targeted specifically within the CGM routine. The background field 
values were taken high enough to achieve convergence to the required accuracy - usually 
values within ±2 sufficed. 

Two different methods were tried to extract the continuum limit. 



Method I 

The first method is a conventional double scaling analysis, as used by Crewther and 
Hamer |T^. Firstly, a sequence of finite-lattice results is extrapolated to the bulk limit, 



^ cxD, at a fixed coupling x (or lattice spacing a). This was done by means of a sequence 
extrapolation routine P2|,p3, namely the alternating VBS algorithm P3]. Errors were esti- 



mated by examining the consistency of the results using different VBS parameters a. The 
sequence displays good convergence to the bulk limit, provided the lattice spacing is not too 
small. An example is given in Table I. 

The estimates of the bulk limit at finite couplings are now extrapolated to the continuum 
limit, a — >^ 0, by performing a polynomial fit to the data in powers of or ga^ which 

conforms to the expected asymptotic series behavior in the weak coupling limit [|12|. Linear, 
quadratic and cubic fits were made to adjacent sets of data points, and extrapolated to 
a = 0. The convergence between these different estimates as the coupling was decreased was 
used to estimate confidence limits for the result. An example of the extrapolation is shown 
in Figure 3 (see later). Note that at small couplings the estimates of the bulk limit become 
unreliable and carry large error bars. 



Method II 

A second method which was tried is a modified version of an approach due to Irving 
and Thomas [|T^. The method involves a simultaneous scaling of the results with lattice 
size A^ and coupling ga, in the fashion of a "phenomenological re-normalization group". 
Eigenvalues are calculated for a sequence of lattice sizes A^ and couplings ga, constrained by 
the requirement 

Nga = c (a constant) (42) 

where c is a constant, so that as A^ — > oo, (7a — > 0, approaching the continuum in a single 
limiting process. The limit of the sequence can again be estimated using the VBS routine. 
For the massless case, m/ g = 0, this procedure works extremely well, as discovered by Irving 



and Thomas JTSl, and shown in Table II. For m/ g = 0, we find by this procedure 
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Mi/g = 0.56417(2) (43) 
Ms/c/ = 1.1284(2) , (44) 



to be compared with the exact results 



Mi/g = 0.564189583 (45) 
M2/g = 1.128379167. (46) 



Irving and Thomas in their previous study with = 14 were able to match and 
even surpass the results obtained with our present Method I. 

For the massive case,however, this procedure turns out to be rather less successful. The 
reason is not hard to find. For a finite system, we have: 

c = Nga = gL (47) 

where L is the physical length of the system. To get the correct result for the excited 
state energies, a box size L will have to be chosen which is larger than the physical size 
of the bound state under consideration. For the solvable massless case, the physical size is 
zero. But for the massive case, this turns out to require such large values of L that the 
extrapolation to the continuum limit becomes less accurate than Method I. Therefore this 
method was dropped from further consideration. 



IV. RESULTS AND ANALYSIS 

Figures 1 and 2 show examples of our estimates of the bulk limit of the vector and scalar 
bound-state energies as a function of lattice coupling ga, at small fermion mass {m/g = 
0.125) and large fermion mass {m/g = 16). It can be seen that the behavior of the vector 
energy is monotonic, while the scalar energy shows a pronounced peak at finite coupling, 
which moves inwards to smaller values of ga as {m/g) is increased. 

At small m/g, the vector energy is almost perfectly linear iny = ga, and an extrapolation 
can be made to the continuum limit y = which is accurate to about 0.5%. At larger masses 
there is more curvature in the data, and any structure moves to smaller values of y, as noted 
above in Figure lb). Figure 3 graphs three different estimates of the continuum limit for 
m/g = 32, as functions of y, obtained by making linear, quadratic and cubic fits to the 
data at that value of y. The different estimates converge nicely down to y = 0.1, where our 
bulk data points begin to lose accuracy, leading to a final estimate of the continuum limit 
accurate to about 1%. 

For the scalar state, our estimates of the bulk energy do not reach far inside the peak 
region before they begin to lose accuracy, as seen in Figures 2a) and 2b). Thus the ex- 
trapolations to the continuum limit are substantially less accurate again, at about the 2-3% 
level. 

Table III lists our estimates of the bound-state energies Ei/ g, E2/ g in the continuum 
limit, together with the earlier finite-lattice estimates of Crewther and Hamer [|r^, the more 
recent light-cone estimates of Mo and Perry [0], and the 'fast-moving frame' estimates 



of Kroger and Scheu |T^. It can be seen that our present estimates are 5-10 times more 
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accurate than the old results of Crewther and Hamer at small m/g, but only 2 times 
more accurate at large m/g\ while for the scalar state they are 2-5 times more accurate. 
The other authors have not attached error bars to their data. The higher-order light cone 
results of Mo and Perry [jl6l are in excellent agreement with ours, over the entire range of 
m/ g. Those of Kroger and Scheu |19| lie a little lower than ours for the vector state at both 



small and large mass ends of the spectrum. Harada et al. IjT^ obtain accurate results for 
low mass using a variational low-mass expansion, but because of the O(m^) errors as the 
results approach m/g = 1 the errors start to escalate. 

Figures 4 and 5 show graphs of our results as functions of {m/g), compared with the 
analytic perturbation theory predictions discussed in Section II. At small m/g, the numer- 
ical results match on well with the small- mass series expansions, equations ([To|) and (|12|), 
although the scalar series begins to diverge away at around m/g ~ 0.1 - presumably further 
terms in the series are required to get a reasonable result. At large m/g, however, a definite 
discrepancy seems to occur: the numerical results lie significantly above the series predic- 
tion, by about three standard deviations. This discrepancy is somewhat puzzling. A closer 
look at Figures lb) and 3), shows a tendency for the data to curve downwards at small y, 
which may indicate that our continuum estimates are slightly too high at these large values 
oim/g. As for the results of Mo and Perry it may be that the fermion self energy term 
has not been included in their analysis. 

The results of Kroger and Scheu [jl9l agree extremely well with the non-relativistic series 
at large m/g. At small m/g they fall below both the analytic small-mass series and the 
numerical results of this work and Mo and Perry |T^. This is presumably because they 
have not included the contributions of multi-quark sectors, qqqq and higher, which have 
significant effect in this region. 

For the scalar state, all results agree well with each other, and with the non-relativistic 
series in the regime of large m/g. 



V. CONCLUSIONS 

In this work a new finite lattice calculation of the lowest lying bound state energies in 
the massive Schwinger model was carried out. Lattice sizes up to 22 sites were treated using 
an equal time Hamiltonian lattice formulation. The results were 2 — 10 times more accurate 



than a previous study |12|. The non-relativistic series were calculated to a higher order. 



These results were compared with the light cone estimates of Mo and Perry |]TB[ and the 



"fast moving frame" estimates of Kroger and Scheu [0. The results generally agree very 
well, showing light cone techniques are powerful and effective for two-dimensional models, 
as first demonstrated by Eller et. al. |15|. The numerical results are generally in satisfactory 



agreement with analytic series results at large and small values of m/g, except for a small 
discrepancy in the non-relativistic region for the vector boson. Possible reasons for this 
discrepancy were discussed. 

Calculations of this sort provide a useful demonstration that both lattice gauge theory 
and light cone techniques can give correct and detailed estimates for the behavior of con- 
tinuum quantum field theory, in a case where analytic results are available for comparison. 
The Schwinger model also exhibits some fascinating physical effects in its own right [Q. 
In future work, we hope to apply both finite lattice and density matrix re-normalization 
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group (DMRG) techniques to a study of chrial symmetry and background field effects on 
this model 
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TABLES 



TABLE I. An example of convergence to the bulk limit, for the vector state energy, Ei/g, at 
fixed coupling, ga = 0.3, and fcrmion mass, m/g = 1. The left hand column gives the finite lattice 
results as a function of lattice size N running from 4 to 22. Subsequent columns give sequence 
extrapolations obtained with the VBS algorithm 



N 

4 -0.691618777 

6 -0.310147879 2.369742591 

8 0.023788597 1.352070697 0.558817207 

10 0.290637884 0.681116152 0.541690825 0.578989183 

12 0.449156794 0.517228742 0.510023826 0.509662896 0.509710838 

14 0.496778711 0.509548577 0.509666964 0.509718226 0.509714839 

16 0.506848384 0.509667696 0.509726825 0.509714344 

18 0.509051005 0.509706664 0.509711057 

20 0.509556263 0.509710588 

22 0.509674480 



TABLE IL Convergence of estimates for the vector state energy, Ei/g, at m/g = using 
the Irving-Thomas method (Method II). The left hand column gives the finite-lattice results as a 
function of N running from 4 to 22. Subsequent columns give sequence extrapolations obtained 
with the VBS algorithm 

N 



4 


0.596129337 










6 


0.581212680 


0.571805670 








8 


0.575443766 


0.569350750 


0.564406265 






10 


0.572480497 


0.568058281 


0.564308303 


0.564213492 




12 


0.570706177 


0.567268711 


0.564260122 


0.564177411 


0.564183130 


14 


0.569535913 


0.566739556 


0.564229677 


0.564184816 


0.564172549 


16 


0.568710910 


0.566361604 


0.564211540 


0.564216832 




18 


0.568100325 


0.566078847 


0.564219012 






20 


0.567631383 


0.565859979 








22 


0.567260599 
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TABLE III. Comparison of bound-state energies Ei/g, E2/g as functions of m/g. The fi- 
nite-lattice estimates obtained in this work are compared with earher finite-lattice estimates of 
Crewther and Hamer [12|, light-cone estimates of Eller et al. [^] and Mo and Perry |16|, and the 



results of Kroger and Scheu ||19| 





Method I 




F.llpr et al. 


Mn & Perry 


Kmger & Scheu 


m/g 


this work 


PI 


N 


M 










Vector state 









0.563(1) 


0.56(1) 








0.125 


0.543(2) 


0.54(1) 


0.58 


0.54 


0.528 


0.25 


0.519(4) 


0.52)1) 


0.53 


0.52 


0.511 


0.5 


0.485(3) 


0.50(1) 


0.49 


0.49 


0.489 


1 


0.448(4) 


0.46(1) 


0.45 


0.45 


0.445 


2 


0.394(5) 


0.413(5) 


0.40 


0.40 


0.394 


4 


0.345(5) 


0.358(5) 


0.34 


0.34 


0.339 


8 


0.295(3) 


0.299(5) 


0.28 


0.29 


0.285 


16 


0.243(2) 


0.245(5) 


0.23 


0.24 


0.235 


32 


0.198(2) 


0.197(5) 


0.20 


0.20 


0.191 








Scalar state 









1.11(3) 


1.12(5) 








0.125 


1.22(2) 


1.11(5) 


1.35 


1.16 


1.314 


0.25 


1.24(3) 


1.12(5) 


1.25 


1.19 


1.279 


0.5 


1.20(3) 


1.15(5) 


1.19 


1.17 


1.227 


1 


1.12(3) 


1.19(5) 


1.13 


1.12 


1.128 


2 


1.00(2) 


1.10(5) 


0.98 


0.99 


0.991 


4 


0.85(2) 


0.93(5) 


0.84 


0.84 


0.837 


8 


0.68(1) 


0.77(5) 


0.69 


0.70 


0.690 
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FIGURES 

FIG. 1. Estimated vector state energies Ei/g in the bulk limit [N oo), as a function of 
coupling ga. a) for m/g = 0.125, b) for m/g = 16. The dashed lines are merely to guide the eye. 

FIG. 2. Estimated scalar state energies in the bulk limit {N — > cxd), as a function of 

coupling ga. a) for m/g = 0.125, b) for m/g = 16. The dashed lines are merely to guide the eye. 

FIG. 3. Convergence to the continuum limit for the vector state energy, Ei/g, at fermion mass 
m/g = 32. Points marked with a circle, square and diamond represent linear, quadratic and cubic 
fits respectively. These fits give extrapolated bulk energies at the axis ga = 0, obtained from 
clusters of two,three or four data points around the indicated values of y = ga. The dashed lines 
give the confidence limits of our final estimate. 



FIG. 4. Continuum estimates of the vector state energy, Ei/g, as a function of log2{m/g). The 
open squares represent this work, the diamonds are results of Mo and Perry [16|, the triangles 
represent Kroger and Scheu The long dashed line, dotted line and dot-dash line represent the 
leading order and higher order non-relativistic series, and an expansion around the massless limit, 
respectively. 



FIG. 5. Continuum estimates of the scalar state energy, E2/g, as a function of log2{m/g). The 
open squares represent this work, the diamonds are results of Mo and Perry [16|, the triangles 
represent Kroger and Scheu [^. The long dashed line, dotted line and dot-dash line represent the 
leading order and higher order non-relativistic series, and an expansion around the massless limit, 
respectively. 
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